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ABSTRACT 
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z ~ 1.2 from the initial data of the DEEP2 Redshift Survey. Assuming a mild 
luminosity evolution, the number of companions per luminous galaxy is found to 
evolve as (l + z) m , with m = 0.51 ±0.28; assuming no evolution, m = 1.60 ±0.29. 
Our results imply that only 9% of present-day L* galaxies have undergone major 
Q-f mergers since z ~ 1.2 and that the average major merger rate is about 4 x 10 _/| 
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h 3 Mpc 3 Gyr 1 for z ~ 0.5 — 1.2. Most previous studies have yielded higher 
values. 

Subject headings: galaxies:interactions - galaxies:evolution - large-scale structure 
of Universe 



1. INTRODUCTION 

Galaxy interactions and mergers are an integral part of our current paradigm of the 
hierarchical formation and evolution of galaxies. Such processes are expected to affect the 
morphologies, gas distributions, and stellar populations of galaxies (e.g., Mihos & Henquist 
1994, 1996; Dubinski, Mihos, & Henquist 1996). Although mergers are rare today (Patton et 
al. 2000, hereafter P2000), cold dark matter iV-body simulations show merger rates of halos 
increasing with redshift as (1 + z) m , with 2.5 < m < 3.5 (Governato et al. 1999; Gottlober, 
Klypin, & Kravtsov 2001). 

Observations, however, have yielded results with < m < 4 (Zepf & Koo 1989; Burkey, 
Keel, & Windhorst 1994; Carlberg et al. 1994; Woods, Fahlman, & Richer 1995; Yee & 
Ellingson 1995; Patton et al. 1997; Neuschaefer et al. 1997; Le Fevre et al. 2000; Carlberg et 
al. 2000; Patton et al. 2002; Bundy et al. 2004). The diverse results are likely due to different 
pair criteria, observational techniques, selection effects, and cosmic variance (Abraham 1999; 
P2000). To identify close pairs, the most secure method is via spectroscopic redshifts for 
both galaxies to find kinematic pairs (Patton et al. 2002, hereafter P2002). This Letter 
adopts the approach of P2002 and uses the early data of the DEEP2 (Deep Extragalactic 
Evolutionary Probe 2) Redshift Survey (Davis et al. 2003; S. M. Faber et al. 2004, in 
preparation) to derive the pair fraction and merger rates up to z = 1.2. In §2, we describe 
the sample and selection functions. In §3, we use both the projected separation and the 
relative velocity to select close galaxy pairs and to determine their evolution. Major merger 
rates out to z ~ 1.2 are computed in §4. The results are discussed in §5. Throughout this 
Letter, we adopt a cosmology of H = 70 km s" 1 Mpc" 1 , VL m = 0.3, and Q,\ = 0.7. 

2. DATA AND SELECTION FUNCTIONS 

The DEEP2 Redshift Survey (DEEP2 for short) will measure redshifts for ~ 50, 000 
galaxies at z ~ 1 (Davis et al. 2003) using DEIMOS (DEep Imaging Multi-Object Spectro- 
graph; Faber et al. 2003) on the 10 m Keck II telescope. The survey convers four fields, 
with Field 1 (Extended Groth Strip) being a strip of 0.25 x 2 deg 2 and Fields 2, 3, and 4 
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each being 0.5 x 2 deg 2 . The photometry is based on BRI images taken with the 12k x 8k 
camera on the Canada- France-Hawaii Telescope (see Coil et al. 2004 for details). Galaxies 
are selected for spectroscopy using a limit of Rab — 24.1 mag. Except in Field 1, a two-color 
cut is also applied to exclude galaxies with redshifts z < 0.75. A 1200 line mm" 1 grating 
(R ~ 5000) is used with a spectral range of 6400 - 9000 A, where the [O II] A3727 doublet 
would be visible at z ~ 0.7 — 1.4. The data are reduced using an IDL pipeline developed at 
UC-Berkeley (J. A. Newman et al. 2004, in preparation). The K-correction is derived from 
spectra of local galaxies (C. N. A. Willmer et al. 2004, in preparation). The data used here 
are from Fields 1 and 4, cover ~ 0.4 deg 2 , and contain ~ 5000 galaxies. 

To measure the spectroscopic selection functions (see L. Lin et al. 2004, in preparation, 
for details), we compared the sample with successful redshifts to galaxies in the full photo- 
metric catalog that satisfy the limiting magnitude and two-color cuts. Following analogous 
approaches by Yee, Ellingson, & Carlberg (1996) and P2002, we calculate the spectroscopic 
weight w of each galaxy as l/S, where S is the spectroscopic selection function derived from 
the R magnitude, B — R and R — I colors, R— band surface brightness fin, and local galaxy 
density of the galaxy itself. To correct for bias due to slit collisions (Davis, Gerke, & New- 
man 2004), we also compute the angular weight, w(8), as a function of angular separation 9 
(P2002). For 9 < 20", we find 0.95 < w(9) < 1.5, a result confirmed with the DEEP2 mock 
catalog of Yan, White, & Coil (2004). 



3. PAIR STATISTICS 

A major limitation on the direct measurement of merger fractions for galaxies is the 
difficulty in identifying on-going mergers, especially for distant galaxies. One alternative is 
to count only the pairs with projected separations Ar and relative line-of-sight heliocentric 
velocities Av less than r max and v max , respectively. A large fraction of pairs with physical 
separations less than 20 h" 1 kpc and velocity difference less than 500 km s" 1 appear to 
have disturbed morphologies or signs of interactions, and these galaxies are expected to 
merge within 0.5 Gyr (P2000). In our work, close pairs are defined such that their projected 
separations satisfy 10 h^ 1 kpc < A r < r max and their rest-frame relative velocities Av 
are less than 500 km s -1 . We adopt an inner cutoff with a projected distance 10 h^ 1 so 
as to avoid the ambiguity between very close pairs and single galaxies with multiple star- 
forming knots. We choose values of r max = 30, 50, and 100 h^ 1 kpc, where 30 h^ 1 kpc is 
most likely to include genuine merger pairs, while the two larger separations provide larger 
samples and thus better statistics. To ensure the selection of the same types of galaxies at 
different redshifts in the presence of luminosity evolution, P2002 adopted a specific range 
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in evolution-corrected absolute magnitude M|, defined as M B + Qz, where the evolution is 
parameterized as M(z) = M(z = 0) — Qz. Following P2002, we adopt Q = 1 as a primary- 
model, and we also study the effect of different models on the pair statistics. We restrict 
the analysis to galaxies with luminosities — 21 < < —19 for z = 0.45 — 1.2. Since there 
is only a 2 mag range in our sample, and assuming a constant M/L ratio, most observed 
pairs are thus major mergers, i.e., mass ratios between 1 : 1 and 6:1. Using data from the 
DEEP2 Fields 1 and 4 helps us reduce the effects of cosmic variance. Following P2000 and 
P2002, we compute the average number of companions per galaxy 

^tot 

where N tot is the total number of galaxies within the chosen magnitude range, Wj is the 
spectroscopic weight for the jth companion belonging to the ith galaxy, and w(6)ij is the 
angular weight for each pair as described in Section 2. Details of the weighting scheme can 
be found in Section 5 of P2002 and references therein. This quantity N c is similar to the 
pair fraction when there are few triplets or higher order iV-tuples in the sample, which is the 
case here. In this work, N c will sometimes simply be referred to as the pair fraction. 

Figure 1 shows N c versus redshift z for a sample with v max = 500 km s _1 and 10 hr x kpc 
< Ar < r max with r max = 30, 50, and 100 hr x kpc, respectively. In the case of using r max = 
50 h" 1 kpc, we find 79 paired galaxies out of 2547 galaxies. The derived N c is ~ 8% at z ~ 0.6 
and increases to 10% at z ~ 1.1. Figure 1 also shows results from the SSRS2 (Southern Sky 
Redshift Survey; P2000) at z ~ 0.015 and CNOC2 (Canadian Network for Observational 
Cosmology; P2002) at z ~ 0.3, after corrections that adopt the same cosmology and the 
same luminosity range. 

The DEEP2 sample is i?-band-selected and hence, for redshifts greater than 0.8, is 
biased toward galaxies that are bright in the rest-frame UV, i.e., against faint red galaxies, 
especially at redshifts z > 1.0. To avoid this bias, we divide our sample into blue and 
red galaxies using the color bimodality feature and then repeat the calculation of N c using 
r m ax = 50 h- 1 kpc for blue galaxies at z = 0.3 (CNOC2), 0.6, and 0.85 bins (DEEP2). When 
using only blue galaxies, the values of the pair fraction N^ lue decrease at z — 0.3 and 0.85 
but increase at z — 0.6 compared to the results using the original sample. The changes in 
the measured value of iV c could be up to a factor of 2, depending on the field and redshift. 
To mitigate the underestimation of the pair fraction in the redshift bin 1.0 < z < 1.2, we 
calculate a corrected N c at z — 1.1 as 

N™ r = m&x[N c ,2N b c lue }, (2) 

shown as open triangles in Figure 1. Parameterizing the evolution of N c as JVc(0)(l + z) m , 
we find that the best fit of (N c (0),m) is (0.029 ± 0.005, 1.08 ± 0.40) with X 2 = 0.27 for 
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rmax = 30 h' 1 kpc, (0.068 ± 0.008,0.51 ± 0.28) with X 2 = 0.89 for r max = 50 h' 1 kpc, and 
(0.177 ± 0.014, 0.47 ± 0.18) with X 2 = 1-2 for r max = 100 h' 1 kpc. For the highest redshift 
bin, N^ or instead of N c is used for the above fitting procedure. The best fits are shown as 
long-dashed curves in Figure 1. For three different choices of r max , we find only a small 
increase of N c from redshifts z = to z = 1.2. In figure 1, we also list the upper limits 
of errors from the cosmic variance for each redshift bin calculated by assuming a spherical 
geometry of the survey volume. Incorporating the cosmic variance into the fitting procedure 
slightly raises m to 0.62 ± 0.49 for r max = 50 h" 1 kpc. 

The adopted choice of Q = 1 in the luminosity evolution model sets the luminosity range 
of the sample at each redshift and thus affects our results. To assess this effect, we repeat 
the pair analysis for other choices of Q. The luminosity ranges are chosen such that they 
are locked to —22 < M B < —20 at z — 1. Using r max = 50 hr x kpc and an upper redshift 
limit of z = 1, we find that m varies from 1.60 ± 0.29 for Q = 0, 0.86 ± 0.29 for Q = 0.5, 
0.41 ± 0.30 for Q = 1, to -0.24 ± 0.35 for Q = 2. In figure 2, we plot N c as a function 
of absolute 5-band magnitude (Mb) for three redshift samples. A clear dependence of N c 
on Mb is evident. The trend that m decreases with Q can be understood as the result of 
including fainter galaxies and thus having larger iV c at lower redshifts when adopting higher 
Q values. Nevertheless, even with the lowest choice of Q = to maximize m, models with 
m > 3.5 are ruled out at a 3 a level of confidence for all three choices of separation limits. 



4. MAJOR MERGER RATES 

The comoving merger rate, usually defined as the number of mergers per unit time per 
comoving volume, can be estimated as 

N mg = 0.5n(z)N c (z)C mg T- l g , (3) 

where T mg is the time-scale for physically associated pairs to merge, C mg denotes the fraction 
of galaxies in close pairs that will merge in T mg , and n(z) is the comoving number density of 
galaxies. The factor 0.5 is to convert the number of galaxies into the number of merger events. 
The merger rate calculated above, however, is not suitable for comparison to the merger rate 
derived from morphological approaches. The reason is that while we have restricted the 
luminosity range of companions to compute N c , the morphological approach does not. To 
correct for this restriction, equation (3) can be modified : 

N mg = (0.5 + G)n(z)N c (z)C mg T m 1 g , (4) 

where the added parameter G accounts for the excess number of companions failing to 
fall into our sample. Assuming that the maximum mass ratio needed to yield significant 
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morphology distortions is 4 : 1, we calculate the value of G as 1.24 for —21 < M% < —19 
using the local luminosity function (Driver & De Propris 2003). T mg depends on the relative 
mass ratio, dynamical orbit, and detailed structure of the two merging galaxies, and thus 
it varies from case to case. Here we adopt T mg = 0.5 Gyr, a value suggested by iV-body 
simulations and simplified models (Mihos 1995; P2000). The value C mg is approximately 0.5 
for close pairs with 5 h^ 1 kpc < Ar < 20 Yr 1 kpc and v max = 500 km s _1 ; this estimate 
is based on morphological studies of local close pairs (P2000). Nevertheless, both T mg and 
C mg remain uncertain. 

The derived merger rates are displayed in Figure 3. Here we have applied the pair 
fraction evolution of m — 0.51 derived for the r max = 50 h -1 kpc case to the z ~ (SSRS2) 
N c result for pairs with 5 h' 1 kpc < Ar < 20 hr l kpc. The parameter n(z) is calculated 
as the sum of the number of galaxies in the adopted magnitude range, each weighted by its 
spectroscopic weight and divided by the comoving volume occupied by the included sources. 
The errors shown for DEEP2 measurements represent 40% variations that are typical for close 
pair counts in our sample, and they do not include the uncertainties of T mg and C rng . Also 
shown are the results from Conselice et al. (2003), who relied on morphologically identified 
mergers for Mb < —19 using Hubble Space Telescope (HST) data in the Hubble Deep Field- 
North (filled triangles). Clearly, the comoving merger rate in our data changes little with 
redshift. From z ~ 0.5 to z ~ 1.2, the average is N mg ~ 4x 10~ 4 h 3 Mpc~ 3 Gyr -1 for 
—21 < M% < —19. This value is about 1 order of magnitude lower than the average merger 
rate for galaxies with M B < —19 derived by Conselice et al. (2003). 

Finally, we calculate the merger remnant fraction, f rem , defined as the fraction of 
present-day galaxies that have undergone major mergers (P2000). Adopting the merger 
fraction C mg to be 0.5 and T mg to be 0.5 Gyr, we estimate, using equation (32) in P2000, 
that about 9% of present L* luminous galaxies have undergone major mergers since z ~ 1.2. 

5. DISCUSSION 

Studies up to z ~ 1 using pair counts have found a wide range in the evolution of merger 
rates. For example, based on roughly 300 Canada-France Redshift Survey galaxies measured 
with HST, Le Fevre et al. (2000) concluded that the pair fraction evolves with m = 2.7±0.6 
while the fraction of merger candidates evolves with m = 3.4 ± 0.6. In contrast, Carlberg 
et al. (2000) measured the mean fractional pair luminosity from z = 0.2 to 1 and found no 
evolution with m ~ ± 1.4, consistent with our result. Since the pair fraction depends on 
adopted luminosity limits, differences among studies may depend on the choice of luminosity 
evolution models. Le Fevre et al. (2000), e.g., applied no corrections while Carlberg et al. 



-7- 



(2000) adopted the Q=l model, as we do here. The key advantages of our pair sample over 
previous surveys include having measurements at z > 1, a larger sample, and more restrictive 
pair criteria than those adopted by Le Fevre et al. (2000) and Carlberg et al. (2000). 

The discrepancy in the derived merger rates at z ~ 1 between our estimates and those 
from Conselice et al. (2003) is more difficult to reconcile. Besides the choice of luminosity 
ranges and the uncertainties in C mg and T mg , other factors may explain the discrepancy. 
First, cosmic variance can always come into play, since most morphological studies have been 
forced to use the few small fields covered deeply by HST imaging, whereas pair-count surveys 
cover much larger areas. Second, the two approaches may sample interacting and merging 
galaxies at different stages, with morphological approaches identifying very advanced mergers 
and merger remnants, while the pairs detect some of the same systems before distortions 
are discernible. Matching the derived merger rates from these two approaches requires 
knowledge of the precise timescales of close pairs and of the duration time for the appearance 
of distorted morphologies. However, neither is well understood yet. Third, while our close 
pairs identify only major mergers, morphological criteria may be sensitive to minor mergers 
as well. Finally, both the pair-count and morphology methods are subject to systematic 
uncertainties that require detail modeling (Bell 2004). Given the current discrepancy in 
results, it behooves us to study the connection between kinematic pairs and morphologically 
disturbed galaxies at various redshifts, both via observational approaches (P2000; J. Lotz 
et al. 2004, in preparation) and through numerical simulations (T. J. Cox et al. 2004, in 
preparation). 

The merger rates from pair counts are also an excellent test for those estimated from 
iV-body simulations, since the former reveal the behavior of luminous baryons while the 
latter reflect the behavior of dark halos. Gottlober et al. (2001) defined merger rates in 
iV-body simulations as equivalent to the major merger fraction per gigayear by tracing the 
formation and evolutionary history of each halo. These fractions will differ from the pair 
fractions roughly by a constant, and therefore their redshift evolution is amenable to direct 
comparisons. They found merger rates evolving as (1 + z) 3 up to z ~ 2. This theoretical 
value of m = 3 is significantly higher than the m ~ 0.5 we find from observed pair counts. We 
should, however, be wary about this comparison for the following reasons. First, definitions 
of merger rates in observations and simulations may not be consistent with each other, since 
the halos and visible galaxies span different size scales and since iV-body simulations also 
suffer limitations in resolution. Second, the merger rates/fractions of halos are likely to be 
a function of halo mass, which is suggested by our finding that the galaxy pair fraction is 
a function of luminosity. Pair count works using i^-band-selected samples (e.g., Bundy et 
al. 2004) can also provide another avenue for testing the mass dependence of merger rates 
since the K-band luminosity is more representative of any underlying stellar mass and suffers 
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less from evolution corrections. Although preliminary, mock catalog simulations by E. Van 
Kampen et al. (2004, in preparation) predict flat slopes for the pair-count fractions (see 
Fig. 3 of Bell 2004), just as seen in our observations. We expect a dramatic improvement 
in our understanding of merger histories via pair counts after more realistic comparisons to 
simulations are possible and especially after a 10-fold increase in sample size when DEEP2 
is complete in 2005. 
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Fig. 1. — Pair fraction, N c , as a function of redshift z. The filled squares represent mea- 
surements from DEEP2 Field 1, and the open squares represent those from Field 4. The 
crosses mark the results from SSRS2 and CNOC2. The open triangles show the corrected 
pair fraction N^ or for the highest redshift bin. The long-dashed curves are best fits by the 
form (1 + z) m , using DEEP2, SSRS2, and CNOC2 data (see text); the dot-dashed curves 
represent evolution with m = 2.5 for reference. From the top to the bottom, r max = 100 
h^ 1 kpc, r max = 50 h~ l kpc, and r max = 30 h~ l kpc. Note that the bottom panel has a 
different vertical scale. The error bars shown in the plot and used for fitting are calculated 
by bootstrapping. The numbers appearing in the bottom panel indicate the upper limits of 
errors (in units of percentage) from the cosmic variance for each redshift bin. Those in the 
parentheses are for open symbols. 
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Fig. 2. — Pair fraction, N c , as a function of absolute 5-band magnitude (Mb) using r max = 50 
h^ 1 kpc for different redshift samples as indicated. For each data point, iV c is calculated 
using galaxies with absolute £>-band magnitude between Mb — 1 and Mb + 1. 
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Fig. 3. — Comoving volume merger rate, N mg , as a function of redshift z. Estimates from 
the DEEP2 fields using equation (4) are marked with symbols as indicated; the data from 
Conselice et al. (2003) are represented by the filled triangles. 



